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SUMMARY 


In  preparation  for  an  experimental  evaluation  of  the  photogrammetry  technique  for 
measuring  vibration  mode  shapes,  a  cantilever  steel  plate  was  designed.  This  report 
describes  the  theoretical  methods  used  to  calculate  the  low  order  natural  frequencies  and 
mode  shapes  of  the  plate,  and  compares  one  of  these  with  the  results  obtained  by  a 
traditional  experimental  procedure. 
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1  INTRODUCTION 


In  preparation  for  an  experimental  evaluation  of  the  photogrammetry 
method  for  mode  shape  measurement,  a  steel  specimen  consisting 
essentially  of  a  cantilevered  plate,  welded  along  one  edge  to  a 
supporting  member,  was  required.  To  provide  a  rigorous  test  of  the 
photogrammetry  technique,  the  cantilevered  plate  was  required  to  be 
about  one  metre  square,  the  modal  amplitude  to  have  a  maximum  value  of  ± 
1 . 5  mm  and  the  natural  frequency  of  a  multi-node  mode  to  be  about  25  Hz . 
This  report  describes  the  theory  used  to  determine  the  plate  dimensions, 
and  gives  a  comparison  of  the  theoretical  mode  shape  with  that  measured 
by  traditional  methods.  In  a  subsequent  report  this  measured  mode  shape 
will  be  further  compared  with  the  one  obtained  by  photogrammetry. 


2  THEORETICAL  ANALYSIS 

Assuming  that  the  supporting  structure  is  rigid,  the  problem 
reduces  to  the  analysis  of  a  steel  plate  built  in  along  one  edge,  with 
the  other  three  edges  free.  Reference  1  provides  simple  formulae  for 
calculating  the  natural  frequencies  and  mode  shapes  of  the  first  five 
modes  of  a  square  plate  supported  in  this  way.  The  modal  amplitude  may 
then  be  calculated  from  the  known  input  force  and  the  assumed  value  for 
the  damping  coefficient  of  the  material.  The  fourth  mode  was  selected 
and,  in  order  to  give  the  desired  natural  frequency,  the  thickness  of 
the  plate  was  set  at  three  millimetres.  One  metre  square  steel  plate 
was  not  available  in  this  thickness  so  the  specimen  was  made  1022  mm  by 
902  mm,  with  the  shorter  side  being  built  in. 

However  since  the  plate  is  now  not  square  the  simplified  formulae 
of  Reference  1  no  longer  strictly  apply  and  the  appropriate  mode  shapes 
and  natural  frequencies  were  calculated  in  the  following  manner. 


2.1  Modal  Functions 

Figure  1  shows  the  coordinate  system  for  the  cantilevered  plate. 
Define  two  non-dimensional  coordinates  £  and  q  as  follows 

\  -  x/a  ,  q  =  y/b  (1) 

where  a  and  2b  are  the  lengths  of  the  free  and  fixed  edge  respectively. 
Following  Reference  1,  the  plate  mode  shape  is  expanded  in  terms  of  beam 
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functions.  These  beam  functions  satisfy  the  beam  boundary  conditions 
corresponding  to  the  free  and  clamped  edges  but  their  product  does  not 
satisfy  the  plate  boundary  conditions  exactly. 

In  the  x  direction  define  the  function 

<pt  (4)  =  au(cosh  0^5  -  cos  oqi;)  +  bli(sinh  -  sin  a^)  (2) 

where  ai(  the  spatial  frequency,  is  the  i'th  root  of  the  non-linear 
equation 


cosa  cosha  =  -1  (3) 

(See  Appendix  A  for  discussion  on  the  solution  of  this  transcendental 
equation.)  The  coefficients  au  and  bu  are  then  given  by 

_ 2(coshoq  +  cosoq  ) _  ^  4  ^ 

11  (l-e’a‘)  [e®1  (l+2sinai)  +  cosoq  -  sinaj 

b  _  _ 2(-sinhai  +  sina, ) _ 

11  ( l-e_ai)  [ea‘  (l+2sinoq)  +  cosoq  -  sinoq) 

The  plate  is  symmetric  about  the  x  axis,  so  the  modal  functions  in 
the  y  direction  are  separated  into  symmetric  and  antisymmetric  groups. 
The  antisymmetric  functions  are  defined  by 

VajW  =  a2 j  sinh  Pjl  +  b2j  sin  Pjn  (6> 

where  the  spatial  frequency  Pj  is  the  j'th  root  of  the  equation 

tan  p  =  tanh  p  (3) 


and 


a2j  =  - -  (8) 

1  cosPj  sinhPj  +  sinPj  coshPj 

b2j  =  - ^£^2 -  (9) 

cosPj  sinhPj  +  sinPj  coshPj 

Equation  7  has  a  zero  root  (j=l)  and  the  corresponding  modal  function  is 
determined  from  l'Hopital's  Rule  to  be 

Vtl(q)  -  q  (10) 
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The  corresponding  symmetric  functions  are 

V  »j  (n)  =  a3j  cosh  Yjn  +  b3j  cos  y  jt|  (11) 

where  Yj  is  the  j'th  root  of 
tan  y  =  -tanh  y 
and 

1 

a3j  - 

2coshy  j 

b3j  = 

2cosy j 

Equation  12  also  has  a  zero  root,  and  the  corresponding  modal  function 
is  then  simply 

Vsilh)  =  1  (15) 

(See  Appendix  A  for  a  discussion  on  the  solution  of  equations  3,  7  and 

12.) 

A  general  mode  of  vibration  may  then  be  represented  by  the 
following  bi-orthogonal  expansion. 

W(5,q)  -  I  "«i  X  j.3  (  Mij  <Pi(5)  V.jdl)  +  Ni3  <Pi<5)Vsj(n)  I  (16) 

where  W(^,t|)  is  a  normal  mode.  Since  the  plate  is  homogeneous  and  of 
uniform  thickness  all  the  normal  modes  of  the  plate  are  either  exactly 
symmetric  or  antisymmetric  about  the  x  (or  %)  axis.  That  is,  for  any 
one  mode  W(^,i))  then  for  all  i,j  either  M,^  =  0  or  =  0.  In  this 
case  the  mode  of  interest  is  symmetric  so  equation  16  may  be  simplified 
to 


(12) 

(13) 

(14) 


W(^,ti)  =  ir.ilj.i  Naj  <Pi(^)V.j(n)  (17) 

The  set  of  modal  functions  used  is  of  course  not  infinite  but  is 
truncated.  In  this  case  equation  17  becomes 

w($,n)  -  IwJ’h  Nij  <Pi(4)v,j(n)  (18) 
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2.2  Plate  Equation  Of  Motion 

As  stated  in  Reference  1,  the  classical  equation  of  motion  for 
plate  dynamics  is 

d2w 

DV4w  +  a  — y  =  0  (19) 

at 


where 


D  = 


Eh 


12  (IV) 


v  is  Poisson's  ratio,  a  is  the  plate  area  density,  V2  is  the  Laplacian 
operator,  E  is  Young's  modulus,  h  is  the  plate  thickness  and  w(£,r|,t)  is 
the  plate  transverse  motion. 


3x  3y 

i  il  .  i  £_ 

~  a2  3^2  b2  an2 

V4  i  i!_  +  _L --i!_ 

~  a4  as4  a2b2 

For  sinusoidal  motion  we  may  write 

w(5,n,t)  =  W($,n)e1,Bt 


+ 


i  il. 

b4  an4 


(20) 


(21) 


(22) 


Substituting  this  into  equation  19  leads  to 
(DV4  -  <d  2a)  W  =  0 


(23) 


which  is  a  single  scalar  equation  with  unknowns  and  o>  .  If  the 
functions  <p  and  y  satisfied  all  the  plate  boundary  conditions  then 
Galerkin's  method  could  be  used  to  transform  this  scalar  equation  into 
an  eigenvalue  equation,  with  the  eigenvalues  giving  to  and  the  elements 
of  each  eigenvector  being  the  corresponding  values  of  N^.  But  the 
products  <py  satisfy  only  the  geometric  (clamped  edge)  boundary 
conditions,  so  the  Rayleigh-Ritz  method  which  gives  solutions  that 
inherently  satisfy  the  natural  (free  edges)  boundary  conditions  is  used. 
This  provides  estimates  of  the  first  p  (  p  s  nm  )  natural  frequencies 
and  mode  shapes,  with  the  lower  order  modes  being  better  estimated  than 
the  higher  (due  to  truncation  effects)  . 
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The  shaker  has  a  low  armature  mass  (230  grams)  and  its  natural 
frequency  is  approximately  that  of  the  mode  of  interest,  so  the  effect 
of  the  shaker  on  the  relevant  natural  frequency  and  mode  shape  is 
minimal.  The  calculated  mode  shape  is  shown  in  terms  of  its  contours  in 
Figure  2. 


3  EXPERIMENTAL  MEASUREMENT 

The  structural  damping  in  the  plate  is  sufficiently  low 
(approximately  0.1%)  for  a  normal  mode  to  be  excited  with  a  single 
shaker,  provided  the  excitation  is  not  at  a  nodal  line.  An  array  of 
measuring  stations  was  mapped  out  on  the  plate  (Figure  3)  which  was  then 
excited  sinusoidally  at  27.2  Hz  to  produce  the  required  mode.  A 
travelling  accelerometer  was  used  to  measure  the  modal  displacements  at 
these  points,  and  the  beam  functions  of  the  previous  section  were  fitted 
by  a  least  squares  process  to  these  displacements  to  produce  the 
contours  shown  in  Figure  4.  Note  that  both  the  symmetric  and 
antisymmetric  functions  were  utilised  to  allow  for  possible  asymmetry 
of  the  plate.  (See  Appendix  B  for  a  more  general  interpolation  scheme.) 


4  DISCUSSION 

The  natural  frequency  of  the  mode  of  interest  was  calculated  to  be 
23.6  Hz  whilst  the  measured  natural  frequency  was  27.2  Hz.  The  free 
shaker  has  a  natural  frequency  of  25  Hz  and  was  attached  to  the  plate  at 
a  point  of  relatively  low  amplitude  in  this  mode,  so  it  is  not 
responsible  for  the  frequency  discrepancy.  Also  it  is  extremely 
difficult  to  build  a  plate  which  is  truely  "built  in"  along  one  edge  as 
there  is  nearly  always  some  motion  at  this  edge.  This  should  have  led 
to  the  theoretical  procedure  overestimating  the  natural  frequency  rather 
than  underestimating  it.  The  plate  was  tested  in  the  horizontal 
position  which  gave  rise  to  a  static  deflection  of  about  66  mm  along  the 
free  edge.  Also  it  was  welded  to  the  mounting  along  the  "built  in"  edge 
which  distorted  the  plate  somewhat  such  that  it  was  no  longer  planar  and 
this,  coupled  with  the  curvature  resulting  from  the  static  deflection, 
gives  rise  to  a  greater  stiffness  than  that  of  a  flat  plate. 

The  natural  frequency  of  this  mode  (as  of  the  other  modes  of  the 
plate)  is  dependent  on  both  the  overall  size  of  the  plate  and  its  aspect 
ratio  (the  ratio  of  a  to  2b  in  Figure  1) .  If,  due  to  non-planar 
distortion  of  the  plate  near  the  clamped  edge,  the  "effective"  length  of 
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the  cantilever  were  reduced,  the  fourth  natural  frequency  would  rise, 
giving  better  agreement  with  experiment.  For  example,  if  "a"  were 
reduced  by  20%,  the  natural  frequency  of  this  mode  would  rise  by  about 
11%. 


The  two  mode  shapes  are  displayed  in  Figures  2  and  4.  The  latter 
figure  shows  that  the  excited  mode  is  not  symmetric  about  the  x  axis 
confirming  the  non-uniformity  of  the  plate.  Despite  this,  the  main 
region  of  difference  between  the  two  sets  of  contours  is  in  the  area  of 
the  clamped  edge.  In  Figure  2  the  nodal  lines  run  from  the  right  hand 
(free)  edge  to  the  left  hand  (clamped)  edge  whereas  in  Figure  4  the 
nodal  lines  converge  before  reaching  the  clamped  edge.  The  region  of 
difference  is  one  of  small  slope  in  this  mode  and  small  perturbations  in 
amplitude  result  in  significant  changes  in  the  contour  pattern.  The 
distortions  of  the  plate  cited  above  may  also  result  in  a  non-uniformity 
of  the  stiffness  distribution  in  this  region.  The  correlation 
coefficient  calculated  between  the  two  modes  at  the  measurement  points 
of  Figure  3  is  0.946  which  shows  reasonable  agreement. 


5  CONCLUSION 

Plate  theory  was  used  to  design  a  specimen  for  use  in  a 
photogrammetry  experiment.  Following  construction  of  the  specimen,  the 
natural  frequency  of  the  fourth  mode  was  measured  to  be  some  15%  higher 
than  that  predicted.  This  may  be  may  be  due  to  the  distortions 
introduced  into  the  plate  during  fabrication  and  the  static  deformation 
of  the  plate  under  gravity.  Although  this  classical  plate  theory  was 
not  adequate  to  predict  accurately  the  dynamic  behaviour  of  the 
distorted  plate,  the  resultant  natural  frequency  and  mode  shape  of  the 
test  specimen  satisfied  the  original  requirements,  i.e.  to  be  suitable 
for  the  photogrammetry  experiment. 
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Figure  2.  Contours  depicting  calculated  mode  shape 
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Figure  3  Location  of  measuring  points,  X,  and  shaker,  O. 


Figure  4.  Contours  depicting  measured  mode  shape 


APPENDIX  A 


The  three  transcendental  equations  (3,  7  and  12)  given  in  Section 
2.1  must  be  solved  for  the  spatial  frequencies  a  ,  p  and  y  .  This  is 
done  by  using  the  Newton-Raphson  method  with  initial  estimates  suitable 
to  ensure  a  rapid  rate  of  convergence. 

A  sufficient  condition  for  the  Newton-Raphson  procedure  to  converge 
to  a  root  of  the  equation 

F (X) =0  ( A— 1 ) 


is  given  by  the  total  stability  theorem.  That  is,  we  may  linearise  the 
equation 


Xn+1  =  Xn  -  nxn)/F'  (Xn) 


about  X„  and  obtain 


sxntl 


8X_ 


F'<xn>5xn  +  ^l.F"  5Xn  +  hot 
F'(Xn)  (F'(Xn))2 


=  F  (Xn)  F"  (Xn)  5Xn 
(F'  (Xn)  ) 2 


+ 


HOT 


(A— 2 ) 


(A-3) 


where  HOT  indicates  higher  order  terms.  Provided  that  HOT  is  bounded, 
then  the  Newton-Raphson  method  will  converge  if 


|  F  (xn)  F"  (xn>  |  <  !  . 

(F'  (Xn)  ) 2 


(A— 4 ) 


A-l 


Furthermore 


the  non-dimensional  function 


P<Xn) 


nx,)F"(xn)  j 
(F'  <Xn)  ) 2 


(A— 5 ) 


measures  the  exponential  rate  of  convergence  of  the  n'th  estimate,  Xn, 
towards  the  the  correct  value.  For  initial  estimates,  X0,  "sufficiently 
close"  to  a  root  of  equation  A-l,  the  rate  of  convergence  after  n 
iterations,  as  indicated  by  p(Xn),  will  be  be  at  least  as  fast  as  that 
indicated  by  p(X0).  That  is, 


p(Xn)  S  p(X0) 


Consider  firstly  equation  3  of  Section  2.1.  This  may  be  written  as 
F (a)  =  cosa  cosha  +1=0  (A-6) 

F'  (a)  =  -sina  cosha  +  cosa  sinh  a 
and  F" (a)  =  -2sin(a)  sinh (a) 

,  ,  | (cosa  cosha  +  1) (-2sina  sinha)  i  ,, 

...  p(a)  =  | - —I  (A-7) 

(-sina  cosha  +  cosa  sinha) 

Rewrite  equation  A-6  as 

cosa  =  -1/cosha 

then  as  a  ,  cosa  =»  0.  Thus  for  initial  estimates  take 

aoi  =  (2i-l)ic/2  i=l,2,3 .  (A-8) 

Substituting  these  into  equation  A-7  gives 

p(a01)  =  0.731 
p(a02)  =  0.0359 

p  (Ooj)  =  0.0016 

For  larger  values  of  i,  p (aoi) ->  4e~a. 


A-2 


Note  that  after  two  iterations  these  three  convergence  factors  become 

p(a2i)  =  0.0087 
p(<x22)  =  0.23  X  10'6 
p(a23)  =  0.7  X  10'12 

This  indicates  that  convergence  is  rapid. 

Similarly  equation  7  may  be  written  as 

G(P)  =  tanp  -  tanhp  =  0  (A-  9 ) 

G'  (P)  =  tan2p  +  tanh2p 

and  G"  (P)  =  2  (tanp  +  tanhp  +  tan3p  -  tanh3p) 

The  convergence  function  is  again 


P(P)  = 


|G(p)  G"(P)  | 
(G'  (P) ) 2 


and  the  initial  estimates  are  given  by 


P0j  =  (4j-3)ii/4  j=2, 3, 4 - 


(A— 10) 


(A-ll) 


Note  that  the  first  root  of  equation  A-9  is  P  =  0,  so  this  is  not 
obtained  by  iteration,  and  the  corresponding  modal  function  is 


val(n)  =  n 

Substituting  the  initial  estimates  (j=2,3,4..)  into  equation  A-10  gives 

p (P02>  =  °-7?  ^  10"3 

p(P03)  =  0.14  X  10'5 

p(P04)  =  0.27  X  10'8 

For  higher  values  of  j,  p(P0j)  -» 2e’2|i,  again  showing  that  convergence 

is  rapid. 

The  symmetric  functions  from  Section  2.1  lead  to  similar 
expressions.  In  this  case  the  initial  estimates  are  given  by 

Y  oj  =  (4j-5)n/4  j=2, 3,4 -  (A-12) 


A-3 


where  again  the  first  root  is  y  =  0,  and  the  corresponding  modal 
function  is 


vsi(n)  =  i 

The  corresponding  convergence  factors  are 

P  (Y02)  =  0.018 
p(y03)  =  0.34  X  10"4 
P (y 04)  =  0.63  x  10'7 

and,  as  before,  p  (y oj)  ->  2e'2y  for  larger  values  of  j,  showing  rapid 
convergence . 


APPENDIX  B 


As  described  in  Section  3,  beam  functions  with  the  correct  boundary 
conditions  were  used  to  interpolate  the  measured  modal  displacements  to 
sufficient  other  locations  on  the  plate  to  produce  the  contours  shown  in 
Figure  4.  A  more  general  interpolation  technique  which  is  not 
constrained  by  particular  edge  conditions  but  which  fits  a  surface  of 
minimum  curvature  to  the  measured  mode  shape  may  be  used. 

The  equation  which  gives  the  displacement  of  an  infinite  plate 
subjected  to  a  number  (n)  of  point  loads  is 

w(x,y)  =  a  +  bx  +  cy  +  r2  log(rl)  (B-l) 

where  r2  =  ((x-xj2  +  (y-y^2)) 

and  f*  is  the  point  load  applied  at  the  point  (x^yj. 

This  same  equation  results  from  the  solution  to  the  problem  of  fitting  a 
surface  of  minimum  curvature  with  the  constraint  that  it  pass  through  n 
points  (x^y^wj  where  (xi,yl)  are  the  coordinates  of  the  i'th 
measurement  point  and  is  the  measured  modal  amplitude  at  this  point . 

Equation  B-l  may  be  evaluated  at  the  n  measurement  points,  leading 


to  n  simultaneous  equations  in  the  n+3  unknowns  a,b, c, flr .  .  . ,  fn.  The 
three  extra  equations  result  from  the  equilibrium  conditions 

Xi-l  ft  =  0  (B-2) 

I".!  fAXi  =  0  ( B— 3 ) 

E"-l  fiYi  -  0  (B-4) 


B-l 


Note  that  these  last  three  equations  are  independent  of  the  measured 
displacements,  w^ 

Ordering  these  equatirns  in  the  sequence  B-2,  B-3,  B-4  and  then  the 
n  equations  resulting  from  B-l,  we  may  express  the  resulting  n+3 
simultaneous  equations  in  matrix  form  as 

U  =  J  V  (B-5) 

where  U  is  the  vector  of  elements  0,0,0,wlf  . .  .  ,wn 

V  is  the  vector  of  elements  a,b,c,  f1# . . . ,  fn, 

and  J  is  a  symmetric  matrix. 

For  non-trivial  solutions  J  may  be  inverted  to  give 

V  =  J'1  U  IB-6) 

Because  the  first  three  elements  of  U  are  zero  the  matrix  J’1  may  be 
partitioned  after  the  first  three  rows  and  columns,  and  the  lower  right 
n  x  n  submatrix  extracted  to  give 

F  =  K  W  IB-') 

where  F  is  the  vector  of  elements  fx,...,fn, 

w  is  the  vector  of  elements  wt,  ...,wn, 

and  K  is  a  symmetric  submatrix  of  J  . 

Thus  equation  B-6  may  be  solved  for  the  coefficients  in  Equation 
B-l  and  the  mode  shape  evaluated  everywhere  on  the  plate.  However  this 
surface  is  constrained  to  pass  exactly  through  the  measured  mode  shape 
displacements  which  are  subject  to  measurement  noise.  Consequently  it 
is  better  to  formulate  the  solution  in  terms  of  a  "relaxed"  least 
squares  problem  subject  to  a  constraint. 

Returning  to  the  analogy  of  the  plate  acting  under  point  loads,  the 
work  done  in  applying  these  loads  is 

work  =  FT  W  ( B— 8 ) 

where  the  superscript  T  indicates  transposition. 

The  work  done  equals  the  increase  in  potential  energy  of  the  deformed 
plate  and  this  in  turn  is  related  to  the  curvature  of  the  plate.  So  to 
minimise  the  curvature  we  may  minimise  the  work.  Thus  the  problem 
becomes : 
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Minimise  FT  W 
subject  to  the  constraint 
[W-W*  ] T  [W-W*]  =  n  a 

where  F  is  the  vector  of  n  coefficients  fir 

W  is  the  vector  of  "noisy"  measured  displacements, 

W  is  the  vector  of  "true"  displacements  (i.e.  without  measurement 

noise) 

and  a2  is  the  variance  of  the  measurement  noise. 


Using  Lagrange  multipliers  this  may  be  written  as 

0.5  FT  W  +0.5  X  [W-W*]T  [W-w'l  =  E  (B-9) 

Substituting  Equation  B-7  into  Equation  B-9  gives 

0.5  WT  K  W  +0.5  X  [W-W*]T  [W-W*]  =  E  (B-10) 

Application  of  the  Variational  Principle  gives 

(K  +  XI]  W  =  X  W*  (B-ll) 

where  I  is  the  unit  matrix. 

Put  e  =  1/X  < B— 1 2 ) 

Then  Equation  B-ll  becomes 

[I  +  eK]  W  =  W*  (B-13) 


It  is  then  required  to  find  the  unique  e  such  that  the  vector  W  obtained 
from  equation  B-13  satisfies  the  equation 

[W-W* ] T  [W-W*]  -  n  o2  =  0  ( B— 1 4 ) 

This  may  be  done  by  iteration  with  the  initial  value  of  e  being 

^initial  =  Trace  (K)  (B-15) 

In  the  present  case  (a  cantilever  plate)  the  edge  conditions  are 
well  stated  and  the  corresponding  beam  functions  have  already  been 
determined  for  the  theoretical  analysis  of  Section  2.  Consequently 
these  beam  functions  are  best  used  for  interpolation  here.  However  in  a 
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more  general  problem  the  necessary  beam  functions  may  not  be  readily 
available  and  the  application  of  the  approach  in  this  appendix  may  be 
useful . 
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